CERN-TH. 7390/94 
TAUP-2191-94 



Large N phase transitions 
and multi-critical behaviour 
in generalized 2D QCD 



B. Rusakov* and S. Yankielowicz** |^ 

*ICTP, High Energy Section, 
P.O.Box 586, 1-34014 Trieste, Italy 

** School of Physics and Astronomy 
Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel-Aviv University, Tel-Aviv 69978, Israel 
and 

Theory Division, CERN 
CII-1211 Geneva 23, Switzerland 



Abstract. 

Using matrix model techniques we investigate the large N limit of generalized 2D Yang- 
Mills theory. The model has a very rich phase structure. It exhibits multi-critical behavior 
and reveals a third order phase transitions at all genera besides torus. This is to be 
contrasted with ordinary 2D Yang-Mills which, at large N, exhibits phase transition only 
for spherical topology. 
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Recently, generalized 2D pure gauge SU(N) and U(N) models (gYM2) have been 
proposed IH and solved at arbitrary N p. The solution generalizes the one obtained for 
ordinary 2D Yang-Mills model (YM2) §-§. 

In particular, the partition function for a closed surface of genus g and area A (coupling 
constant absorbed into area) has the form: 

ZM)=i:dT''^^^v{-4^Mr)) (1) 

Here, r is an irreducible representation of the gauge group, which is parametrized by its 
highest weight components, ni > ... > Ur^, associated with the lengths of lines in the 
Young table; dr is its dimension, 

dr= n (1+^^) (2) 

l<j<j<Af 

and 



J - « 



Mr) = Nj:^C,ir) (3) 



fc>0 

with Cfc(r) being the k-th Casimir operator eigenvalue. 



N 



Ck{r) = Y.lhr, (4) 



where 



7. = n (1 - k = n, + N-i. (5) 

Formula (|I|) generalizes the known result for YM2 which corresponds to A(r) 

with only the k = 2 term present. The eigenvalue of the second Casimir operator is 

N 

C2{r) = Y^n,{n, + N + l-2i) . (6) 

i=l 

A generalization of the result for Wilson loop averages is, 

m A 

W,{C) = E '^n-r^ n exp ( - ^A(r.)) , (7) 

ri,...,rm k=l 

where m is the number of windows, A^s are the corresponding areas, gk is the "genus per 
window" and the coefficients $r-i...rm are the U(N) (SU(N)) group factors which depend 
on the contour topology (see Ref.[Q for details). The result depicted in (0) is obtained 
from the YM2 result by the replacement C2(r) — ^ A(r) 

In Ref. the large N expansion of the generalized action was carried out and a stringy 
description was derived. The partition function (Q) can be written as a sum of 2D string 
maps, where maps with branch points of degree higher than one as well as "microscopic 
surfaces" play an important role. This generalizes the string interpretation of YM2 given 
previously by Gross and Taylor 0-0. It is still a challenging problem to find the string 
action which gives rise to the sum of maps that reproduces the gauge theory partition 



function. Recently, an important progress toward this goal has been made P],[10 . The 
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results of Ref. [^] reveal the topological aspects of the underlying string theory associated 
with holomorphic maps. 



In the present letter we continue to study the gYM2 theory. We employ the large N 



approach, elaborated in details in |TT|, which is well known within the context of matrix 
models. We demonstrate that gYM2 exhibits a much richer structure than YM2. In 
particular, besides the phase transition realized |T^ in YM2 on a sphere, there are phase 
transitions at any genus in gYM2. Moreover, at any genus the model accommodates 
multi-critical behaviour. 

Following |]lT| we write the partition function (|I|) at large N as a path integral over 
continuous Young tables. Namely, we introduce the continuous function. 



^ ^ N^oc AT V 2 
Then, the partition function takes the form: 

Z = 



X 



k 
N 



Y[ dh{x)e 

0<x<l 



S[h{x)] 



S = ^ £dx (^-AV[h{x)] + (2 - 2g) £ log \h{x) - h{y) 



(8) 

(9) 

(10) 



where 



k>0 



It is implied that the path integral (^) should be taken over h{x)^s which satisfy the 
constraint 

> 1 
dx ~ 

which is reminiscence of the dominance condition for group representations. 
Then, the saddle point equation is 



(^~^ is the surface density of Euler characteristics). 
Introducing the density 

Pih) = $ 



A 



2^7 



dh 



which should be positive and normalized to 



we rewrite (1TT|) as 



dh p{h) 



dxp{x) 



h ~ X 

The constraint h'(x) > 1 now takes the form 



a<h<b 



12) 



(13) 



Pih) < 1 . 



(14) 
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The solution of (|T3|) is 

p{h) = ^P{h)^{h-a){b-h) , (15) 



where the polynomial P{h) and the values a and b are uniquely defined by the h oo 
asymptotic behavior 

^CV'ih) - P{h)^{h-a){h-b) ~ i . (16) 
In the YM2 case {V = /i^), the solution (constraint (0) is ignored) is: 

p{h) = ^Vb^ - \h\<l e&' = 2 . (17) 

TT 



Thus, the solution exists only for genus zero, i.e., the sphere. At 2^ < vr^ Eq.(|14[) is 
satisfied for all h. It means that in sum over representations we can ignore the dominance 
condition and all n^'s now run independently from —00 to +00. In other words the U(N) 
(as well as SU(N)) model now becomes (effectively) 0^U(1) model, i.e., abelian. In terms 
of the original U(N) matrix variables it means that all unitary matrices are frozen near 
I. This is the weak coupling (small area) phase of the model. 

For 2^ > TT^, there is a region of values of h where solution (|T^ does not satisfy 
the inequality (p!4D. In this region the inequality (|1^) is saturated, p{h) = 1. Then, the 
solution of the saddle point equation is 



p{h) = I I ioT -b < h <b 



2 

nah 



a2-/i2)(/i2_52)ni(^-|i,j) ioTb<h<a 



where Ili{x, k) is the complete elliptic integral of the third kind with modulus k = - and 
the parameters are determined by the equations 

a{2E - k''^K) = 1 aA = AK , (19) 

whereQ k' = \J\ — k'^ is the complementary modulus and 

K = Kik) = C '^^ E = Eik) = 

Jo y'(l-x2)(l-Pa;2) 




At the critical value, 2^c = {A^ = vr ), the third order phase transition takes place 



It was observed by Douglas and Kazakov [12 



We observe that in the generalized case the constraint (|T^ can be satisfied also for 
negative ^ by an appropriate choice of the coupling constants a^, which implies the exis- 
tence of phase transitions at higher genera, except the torus. (In the latter case the saddle 
point equation is an algebraic equation for h with solution h = const and, therefore, the 
constraint /i' > 1 is not satisfied). Moreover, there is a multi-critical behaviour at each 
genus. Indeed, since the semicircle distribution ( p!7| ) is now deformed by a polynomial 
(|1|), there are several disconnected regions of h where the constraint {\L^) is satisfied. 

^AU notations are taken from the book lOi. 
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Consider, for example, the simplest generalized case 

v{h) = h' + 



which is the well-known if model. The solution is [15 



p{h) = i[l + ^(a + 6) + ah]J{h-a){b-h) , (20) 

IT 2 * 



and the values of a and b are determined from 

f(6-a)2 + (a + 6)[l + f(a + 6)] = 

Cib ~ aY[l + a{a + b)] =8 , 



(21) 



which leads to 



2a y 5 



h = ^''-^ + 2 



(22) 



with 



o 1 , , 

- -z+ — = . (23) 



4 2^ 

At negative ^ we should take the z < solution which exists for 



2^ e 

a < 



3V12 

The critical value = Cc(a) should be defined from (pO|) through the condition 

Pmax(/i) = 1. (24) 

It is clear that similar situation prevails for the general case. For ^ < ^c, where is 
defined by (p^ and p defined by (p!5[) , the weak coupling phase is realized. 

In the strong coupling phase (^ > ^c) the saddle point equation (|T^) should be solved 
in presence of the constraint (0) . The only solution is 

{1 o,k<h<bk k = l,...,m 

(25) 
p{h) bk<h<ak+i k = 0,...,m 

where feo = Om+i = b and m is the number of intervals (a^, bk) where 

1 



-P{h)J{h-a){b-h) > 1 
7r * 

The function p{h) has to be defined via the equation 



iem-EiogF^ = Er"^ ■ (26) 

2^ ^ ^ fctl /i-flfc t'o^b, h-x ^ ' 
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Wilson loop averages can be calculated following the same steps as in the YM2 case 
on a sphere |jl6|],[|l^. A simple (without self- intersections) contractible loop divides 
the surface into two "windows" with areas Ai,A2 and genera per "window" gi,g2- The 
corresponding expectation value can be calculated in a similar way to what was done for 



YM2 on a sphere []I6 



where the quantities ri{h) and D{h) are precisely the same as in the spherical case (see 



Ref. ||16||) and are expressed through p{h). We note that Eg. ([27|) is symmetric with re- 
spect to the exchange {Ai,gi) (7^2, (72)- Again, as in the spherical case, VF(Ai,0) = 
W{Q,A2) = 1. 

Expression (|27D is valid in both the weak and the strong coupling phases. A complete 
information concerning the phase transition is contained in the function p{h). This func- 
tion is given by (|15[) in the weak coupling phase and by ([25|)-(p6|) in the strong couphng 
phase. 



Starting from ([27| ) the complete set of Wilson loop averages can be calculated by 



means of the loop equation [T^ in its simple two-dimensional formulation |jT9|, similarly 



to what was done for YM2 on a sphere |]T7| 



Several remarks are in order. First, the non-trivial saddle point for fi' > 1 exists only if 
the coupling constant corresponding to the highest Casimir eigenvalue in the generalized 
potential (^ is negative. In terms of the matrix model it implies that the potential 
V{h) has a negative leading term, i.e., we deal with a matrix model with an upside-down 
potential. It is easy to check that there is no phase transition on a sphere for the case 
of an upside-down potential. Our second remark concerns the function = ^c{g)- For 
any particular model (fixed couplings Uk and maximal order of the Casimir operators 
entering in (^) the critical value is fixed (it is unambiguously defined by the saddle 
point equation). Therefore, increasing the genus the critical area increases proportionally 
in order to keep fixed the critical number of handles per area (the density of the Euler 
characteristics). Also, if we fix the area A of the surface, then with the increase of the 
genus the number of multi-critical points decreases and for g > 1 + there is no phase 
transition. 

There are several open questions which follow from our observations. We have pre- 
sented only a qualitative picture of the critical behavior in gYM2. It is of interest to get 
the critical exponents associated with the phase transitions. 



In the weak coupling phase, where the constraint (|l^ can be ignored, the model coin- 
cides with the usual (large N) hermitian matrix model with arbitrary potential. The latter 
has been investigated thoroughly during recent years. It describes (upon fine tuning of 
the parameters) two-dimensional quantum gravity (coupled to non-unitary CFT matter). 
It can be also viewed as string embedded in D = space. It is very intriguing to find out 
whether there exists an hermitian matrix model which corresponds to the strong coupling 
phase of gYM2. If such a matrix model description is found, then the obvious question 
is to which quantum gravity string theory it corresponds and what is its relation to the 
stringy formulation of gYM2 Clearly, the stringy approach to gYM2 holds 

■^An expectation value of non-contractible loop might be zero due to pure topological reason. A 
complete set of zero- valued loops was found in [|| . 
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in the strong coupling phase. As a first step toward the construction of such a matrix 
model one could try to introduce the constraint ([l^) directly into the action (p!0[) of the 
model (the constraint can be included in the path integral over h{x) as a step function 
6(\h'{x)\ — 1) which can be incorporated into the action, for example, via a Lagrange 
multiplier or through some convenient parametrization) . 

One can start with the lattice formulation of gYM2 and integrate systematically over 
the link gauge variables. This gives rise to an effective action for the auxilliary B field 
which is again a matrix model. The matrix model formulation can help toward finding 
an integrable structure of gYM2. 
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